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Periodic  structures  in  the  equatorial  ionosphere 

Cheryl  Y.  Huang,1  Susan  H.  Delay,2  Patrick  A.  Roddy,1  and  Eric  K.  Sutton1 
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[i]  The  plasma  densities  detected  by  the  Planar  Langmuir  Probe  (PLP)  on  board  the 
Communication/Navigation  Outage  Forecasting  System  (C/NOFS)  satellite  have  been 
analyzed  for  the  June  2008  period.  This  interval,  which  corresponds  to  one  of  the 
quietest  periods  in  the  space  era,  exhibited  broad  plasma  decreases  (BPDs)  which  we  have 
reported  previously.  In  order  to  treat  the  data  quantitatively,  we  have  detrended  the  PLP 
data  by  using  the  International  Reference  Ionosphere  model  to  remove  variations  in 
density  due  to  changes  in  spacecraft  altitude  and  latitude  along  the  orbit.  In  this  paper  we 
present  results  of  a  statistical  analysis  of  C/NOFS  detrended  plasma  densities  during 
June  2008  as  well  as  neutral  densities  from  the  Gravity  Recovery  and  Climate  Experiment 
measured  during  the  same  period.  The  results  show  periodic  structures  in  both  charged  and 
neutral  species,  most  evident  in  limited  local  times.  These  structures  resemble  wave-4 
nonmigrating  tides.  BPDs  persist  in  the  detrended  data  and  appear  as  one  of  the  minima 
in  the  periodic  signatures,  most  strikingly  on  the  nightside. 

Citation:  Huang,  C.  Y.,  S.  H.  Delay,  P.  A.  Roddy,  and  E.  K.  Sutton  (2011),  Periodic  structures  in  the  equatorial  ionosphere, 
Radio  Sci.,  46,  RS0D14,  doi:  10.1 029/20 10RS004569. 


1.  Introduction 

[2]  It  is  generally  accepted  that  the  equatorial  low-latitude 
ionosphere  is  influenced  by  solar  radiation,  magnetic  activity 
and  variations  in  solar  cycle  [de  La  Beaujardiere  etal,  2004] . 
These  combine  to  create  a  wealth  of  plasma  interactions 
which  reflect  changes  in  local  time,  latitude,  altitude  and  sea¬ 
son.  Extensive  studies  of  plasma  irregularities  have  estab¬ 
lished  the  complex  relationship  between  E  and  F  region 
electric  fields  and  neutral  winds  driving  the  growth  of  the 
Rayleigh-Taylor  instability  [see  de  La  Beaujardiere  et  al., 
2004,  and  references  therein].  The  C/NOFS  satellite  was 
launched  in  April  2008  with  the  aim  of  predicting  commu¬ 
nications  outages  due  to  these  equatorial  irregularity  structures. 

[3]  However  more  recent  studies  have  indicated  the  impor¬ 
tance  of  the  troposphere  as  a  contributor  to  ionospheric 
structure.  Migrating  and  nonmigrating  tides  [ Chapman  and 
Lindzen,  1970]  are  caused  by  an  interaction  between  tropo¬ 
spheric  latent  heat  and  solar  diurnal  tides  which  are  generated 
by  solar  radiation  [Hagan  and  Forbes,  2002,  2003].  Dis¬ 
covery  that  these  tides  penetrate  deep  into  the  upper  F  region 
was  unexpected,  as  planetary  waves  are  expected  to  damp 
out  far  below  the  F  region  peak.  Results  from  the  C/NOFS 
mission  have  confirmed  that  tidal  structures  are  dominant  in 
the  ionosphere,  being  reported  from  a  number  of  separate 
investigations  over  the  entire  altitude  range  of  the  satellite 
which  extends  to  800  km  [Pfaff  et  al.,  2010;  Araujo -Pradere 
et  al.,  2010;  Dao  et  al.,  2010].  In  this  paper  we  present  an 
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analysis  of  periodic  structures  resembling  tides  found  in 
C/NOFS  plasma  densities,  and  focus  on  one  particular  phe¬ 
nomenon  called  a  Broad  Plasma  Decrease  (BPD). 

[4]  In  our  previous  study  [Huang  et  al.,  2009],  we  reported 
BPDs  observed  on  C/NOFS  as  well  as  on  the  CHAllenging 
Minisatellite  Payload  (CHAMP)  and  Defense  Meteorological 
Satellite  Program  (DMSP)  satellites.  An  example  of  a  BPD 
on  C/NOFS  is  shown  in  Figure  1  which  occurred  on  5  June 
2008,  orbit  734.  Figure  1  (top)  shows  the  plasma  density 
(1  min  averages  in  black,  1  s  averages  in  green)  together  with 
the  satellite  altitude  (in  blue,  with  heavier  trace  when  the 
spacecraft  is  on  the  nightside).  In  this  paper  we  refer  only  to 
the  1  min  averages.  Ephemeris  information  is  given  in  the 
table  at  bottom.  Figure  1  (middle)  gives  an  indication  of  the 
frequency  of  the  fluctuations,  not  referenced  in  this  study. 
Figure  1  (bottom)  shows  the  orbit  in  red,  with  a  heavier  trace 
showing  when  the  satellite  is  on  the  nightside.  Also  shown 
in  Figure  1  (bottom)  is  the  magnetic  equator  in  green.  The 
large  decrease  in  plasma  density  commencing  at  approxi¬ 
mately  04:09  UT  is  the  signature  of  the  BPD.  As  noted  in 
our  past  study,  BPDs  are  a  regular  feature  of  solar  minima, 
with  the  largest  BPDs  seen  during  the  recent  past  solar 
minimum.  They  are  most  evident  at  South  Atlantic  Anomaly 
longitudes,  although  they  are  also  detected  in  other  longi¬ 
tude  regions.  In  general  BPDs  are  detected  between  May 
and  July  on  DMSP  in  every  solar  minimum  year  in  the 
2006-2010  interval  in  the  recent  solar  minimum,  but  also 
in  1996  in  the  previous  solar  minimum  [Huang  et  al., 
2009].  In  order  to  show  BPDs  together  with  periodic 
structures  clearly,  we  restrict  analysis  in  this  paper  to  the  June 
2008  period,  although  they  also  occur  in  June  2009-2010. 

[5]  Quantitative  analysis  of  these  density  variations  is 
made  difficult  by  the  changes  in  C/NOFS  location.  Note  that 
from  the  start  of  the  density  decrease  shown  in  Figure  1  to 
the  end,  the  satellite  altitude  changes  from  approximately 
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Figure  1.  C/NOFS  observations  of  a  Broad  Plasma  Depletion  (BPD)  on  5  June  2008.  (top)  The  plasma 
density  (black  line)  in  1  min  resolution  as  well  as  1  s  resolution  (green  line).  C/NOFS  altitude  is  indicated 
by  the  blue  line,  with  the  heavy  dotted  section  showing  when  the  satellite  is  in  darkness,  (middle)  The 
spectrum  of  the  high  time  resolution  density  fluctuations,  (bottom)  The  magnetic  equator  is  shown  by 
a  green  trace,  and  satellite  orbit  is  indicated  by  the  red  trace,  with  the  heavy  dotted  section  indicating  when 
the  spacecraft  is  in  darkness.  The  BPD  starts  at  04:09  UT,  with  recovery  after  04:33  UT.  Ephemeris  infor¬ 
mation  is  given  along  the  bottom. 


400  km  at  04:09  UT  to  700  km  at  04:33  UT.  In  order  to 
reduce  the  density  variations  due  to  change  in  satellite 
altitude  and  latitude  we  have  detrended  the  plasma  data  by 
using  the  IRI  model.  Altitude  is  the  main  contributor  to 
orbital  variations,  hence  we  refer  to  our  detrending  process 
as  altitude-detrending  for  clarity.  The  data  are  all  referenced 
to  the  equator  at  an  altitude  of  500  km.  We  then  analyze  the 
altitude-detrended  data  in  order  to  quantify  the  density  varia¬ 
tions.  The  focus  in  this  paper  is  on  the  June  2008  period. 
This  period  is  marked  by  extremely  quiet  magnetic  condi¬ 
tions,  with  solar  UV  fluxes  the  lowest  seen  for  many  dec¬ 
ades.  There  are  high-speed  streams  in  the  solar  wind  which 
affect  the  ionosphere  and  thermosphere  periodically,  but  not 
daily.  BPDs  are  observed  both  in  the  presence  and  absence 
of  fast  streams.  Without  external  driving  by  the  solar  wind, 
the  ionosphere  is  affected  strongly  by  internal  processes. 

[6]  We  carried  out  a  statistical  analysis  on  June  2008 
accelerometer  data  from  the  GRACE  satellite  to  determine 


if  the  structures  that  appeared  in  the  plasma  analysis  were 
also  present  at  the  same  time  in  the  neutral  densities.  This 
appears  to  be  the  case,  although  the  coupling  between  ions 
and  neutrals  is  not  straightforward. 

2.  Instrumentation 

[7]  The  PLP  on  C/NOFS  measures  plasma  densities,  elec¬ 
tron  temperatures  and  density  fluctuations  at  rates  ranging  from 
32  to  1024  samples  s  '.We  have  used  1  min  averaged  data  for 
our  analysis.  In  this  paper  we  focus  on  electron  densities  on 
the  nightside.  During  June  2008,  the  period  of  this  study, 
density  depletions  were  observed  over  the  full  range  of  the 
C/NOFS  orbit,  from  400  to  840  km. 

[8]  Neutral  densities  on  the  Gravity  Recovery  and  Cli¬ 
mate  Experiment  (GRACE)  satellite  were  measured  using 
superSTAR  accelerometers,  using  the  techniques  of  Sutton 
et  al.  [2005]  and  Sutton  [2009].  These  measurements  are 
similar  to  those  made  by  the  STAR  accelerometer  on 
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PLP  and  IRI  Densities  along  CVNOFS  Track  Altitude-Detrended  PL P Density  and  ERI  at  500  lan 


Figure  2.  (a)  The  plasma  density  measured  by  the  Planar  Langmuir  Probe  (PLP)  along  the  satellite  track 
on  5  June  2008,  orbits  733-734.  The  results  are  shown  as  functions  of  geographic  longitude,  altitude, 
and  local  time.  Solid  blue  circles  denote  the  PLP  measurements.  The  values  of  the  International  Reference 
Ionosphere  (IRI)  model  for  the  satellite  location  are  indicated  by  red  squares,  (b)  The  result  of  altitude¬ 
detrending  the  PLP  plasma  density  using  the  IRI  model.  Altitude-detrended  and  IRI  plasma  densities  are 
shown  at  an  altitude  of  500  km  and  at  the  equator. 


CHAMP  [Bruinsma  et  al.,  2004;  Sutton,  2009],  with  a  pre¬ 
cision  an  order  of  magnitude  greater  than  the  CHAMP 
instrument.  In  this  study,  5  s  data  were  used.  During  the 
period  of  interest,  the  satellite  is  in  a  443  x  485  km  orbit. 
Factoring  in  the  oblateness  of  the  Earth,  the  minimum  alti¬ 
tude  in  the  World  Geodetic  System  84  (WGS84)  was  460  km 
and  the  maximum  was  506  km. 

[9]  GRACE  makes  30-32  equatorial  passes  per  day;  15-16 
on  the  ascending  portion  of  its  orbit  and  15-16  on  the  des¬ 
cending.  Given  the  rotation  of  the  Earth,  each  of  the  passes 
occurs  over  a  different  longitude  such  that  after  a  day,  there 
are  15-16  data  points  equally  spaced  in  longitude,  per 
hemisphere.  The  repeat  ground  track  period  is  7  days,  after 
which  time  there  are  about  112  data  points  equally  spaced  in 
longitude,  per  hemisphere.  The  GRACE  plots  in  this  paper 
use  28  days,  so  there  are  4  repeat  ground  tracks  in  lon¬ 
gitude,  with  limited  coverage  of  local  times.  In  this  report 
we  include  only  those  local  times  for  which  we  have  suf¬ 
ficient  numbers  of  observations.  We  restrict  the  data  to  the 
same  latitude  range  as  that  covered  by  the  C/NOFS  satellite 
during  June  2008. 


3.  Observations 

[10]  An  example  of  the  observations  that  stimulated  this 
investigation  is  shown  in  Figure  1.  As  already  noted,  the 
BPD  which  spans  16  min  from  04:09  UT  to  04:25  UT  is 
detected  over  an  altitude  range  from  750  to  600  km.  This 
spatial  variation  makes  it  difficult  to  evaluate  the  true  reduc¬ 
tion  in  density.  In  order  to  quantify  the  density  decreases,  we 
use  the  2007  version  of  the  IRI  model  [Bilitza  and  Reinisch, 
2008]  using  the  default  topside  density  model.  PLP  densities 
are  then  processed  as  follows: 

NDet  =  N(PLP)  *  N(IRI  at  500  km)/N(IRI  at  C/NOFS  location), 


where  N(PLP)  is  the  density  measured  by  PLP  along  the 
C/NOFS  orbit;  N(IRI  at  500km)  is  the  IRI  value  at  the 
equator  at  500  km  altitude;  and  N(1RI  at  C/NOFS  location) 
is  the  IRI  value  at  the  same  location  and  local  time  as  the 
PLP  instrument. 

[11]  By  taking  the  ratio  of  two  IRI  model  values,  we  min¬ 
imize  the  intrusion  of  the  model  in  our  altitude-detrending 
scheme.  We  treat  all  the  altitude-detrended  values  as  if 
measured  at  the  equator  at  a  constant  height  of  500  km.  It 
should  be  noted  that  the  very  low  ionospheric  densities 
during  June  2008  are  not  represented  well  by  IRI.  However, 
our  detrending  scheme  does  not  rely  on  absolute  values  of 
density  in  the  model,  only  on  the  relative  densities  at  dif¬ 
ferent  altitudes. 

[12]  Figure  2  illustrates  the  altitude-detrended  densities 
for  parts  of  orbits  733  and  734  covering  most  of  the  interval 
shown  in  Figure  1.  The  plots  show  densities  as  functions  of 
geographic  longitude,  altitude,  and  local  time.  Figure  2a 
shows  the  original  densities  along  the  C/NOFS  orbit  track, 
together  with  the  IRI  model  values  at  the  satellite  location. 
Figure  2b  shows  the  altitude-detrended  densities  together 
with  the  IRI  model  values  at  500  km  at  the  equator.  Note 
that  the  scales  on  the  two  plots  are  different.  Mapping  PLP 
densities  to  500  km  extends  the  density  range  since  the  orig¬ 
inal  orbit  was  close  to  apogee  on  the  dayside.  In  Figure  2b 
the  BPD  can  be  seen  between  -40°W  (or  320°E)  longitude 
and  about  20°E  longitude.  Note  that  there  are  other  depletions 
evident  in  Figures  2a  and  2b  between  -90°W  and  — 150°W 
longitude.  As  we  will  show  below,  periodic  modulations  are 
strikingly  in  evidence  during  June,  with  the  BPD  the  largest 
feature  within  the  periodic  structures. 

[13]  In  the  original  data,  the  density  decreased  from  8  x 
104  cnT3  at  04:06  UT  (near  satellite  perigee  at  400  km) 
to  1.8  x  103  cnT3  at  04:16  UT  near  our  nominal  altitude 
of  500  km.  In  the  altitude-detrended  data,  the  densities 
are  1.8  x  104  cnT3  and  1.7  x  103  cnT3  at  the  corresponding 
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Figure  3.  Altitude-detrended  PLP  plasma  densities  are  shown  as  functions  of  longitude  for  the  day- 
side,  divided  into  2  h  local  time  averages.  The  data  are  averaged  over  3°  wide  longitudinal  bins. 
For  comparison  the  IRI  values  are  included.  Wave-4  periodicity  can  be  seen  in  the  plots  for  08:00- 
10:00  LT,  10:00-12:00  LT,  and  12:00-14:00  LT.  At  dawn  a  two-wave  structure  can  be  observed. 


times.  The  altitude-detrended  density  shows  a  decrease  of 
91%  within  the  BPD  while  the  original  density  decreases 
by  98%.  Altitude  variation  obscures  the  “true”  density 
change.  In  addition,  altitude-detrending  reveals  density 
variations  which  are  not  as  easy  to  identify  as  the  BPD  in 
Figure  1.  These  will  become  apparent  in  our  statistical  anal¬ 
ysis.  Also  evident  in  Figure  2  are  the  differences  between  the 
measured  plasma  density  during  this  extremely  deep  solar 
minimum  and  the  IRI  model  values.  This  is  not  altogether 
unexpected  given  that  solar  flux  values  during  2008  are  far 
below  past  solar  minimum  observations. 

[14]  Figures  3  and  4  show  the  average  altitude-detrended 
densities,  NDet  on  the  dayside  and  nightside,  respectively, 
together  with  IRI  model  values  for  the  month  of  June  2008. 
The  densities  have  been  sorted  into  2  h  local  time  bins  and 
averaged  over  3°  of  geographic  longitude.  The  averages  are 
shown  as  functions  of  longitude  for  each  local  time.  Figure  3 
(left)  shows  the  results  for  06:00-08:00  LT,  08:00-10:00  LT, 
and  10:00-12:00  LT.  Figure  3  (right)  shows  results  for 
12:00-14:00  LT,  14:00-16:00  LT,  and  16:00-18:00  LT.  This 
is  repeated  in  Figure  4  for  local  times  on  the  nightside,  as 
indicated  in  each  panel.  Although  not  shown  in  Figure  3  or 
Figure  4,  the  standard  deviations  of  NDet  within  each  lon¬ 
gitude  bin  can  be  large,  at  times  comparable  with  the  values 
themselves. 

[15]  The  neutral  densities  derived  from  the  GRACE  accel¬ 
erometer  were  analyzed  to  remove  altitude  and  local  time 
variations  as  follows.  During  an  orbit,  the  near-polar  incli¬ 
nation  of  the  GRACE  satellite  allows  for  near-global  latitu¬ 
dinal  sampling.  Each  latitude  is  sampled  twice  during  an 


orbit,  at  two  different  local  times  that  are  approximately 
12  h  apart.  There  are  15-16  satellite  orbits  per  day,  which 
provides  full  longitudinal  coverage  on  the  ascending  and 
descending  portions  of  the  orbit.  Therefore,  an  effective  zonal 
mean  density,  pzm,  can  be  computed  at  each  latitude  and  alti¬ 
tude  for  each  day.  This  zonal  mean  can  then  be  removed  from 
the  density  data  to  yield  a  percent  residual:  1 00*(  p  -  p^/p^. 
These  residuals  are  binned  according  to  latitude  and  longi¬ 
tude,  for  both  ascending  and  descending  portions  of  the  orbit. 

[16]  The  resulting  plots  show  the  longitudinal  structure  of 
the  thermosphere  as  viewed  from  a  satellite  in  a  nearly  fixed 
local  time  frame.  Longitudinal  structures  of  a  certain  wave 
number,  k,  should  be  interpreted  as  a  sum  of  all  of  the  tidal 
components  with  Is  -  nl  =  k,  where  2/rn  is  the  frequency  of 
oscillation  (in  days)  and  s  is  the  zonal  wave  number  from 
classical  tidal  theory  [see  Lindzen  and  Chapman,  1969]. 

[17]  Two  additional  caveats  arise  from  interpreting  these 
plots.  First,  by  removing  the  daily  zonal  mean  density,  the 
residuals  alias  with  the  migrating  tides  such  that  the  resi¬ 
duals  will  be  smaller  when  the  satellite’s  orbital  plane  is  in 
phase  with  the  migrating  tides,  and  larger  when  out  of 
phase.  In  an  effort  to  minimize  this  problem,  we  use  both  the 
ascending  and  descending  density  values  to  compute  the 
daily  zonal  mean  density  for  a  set  of  data  points  at  a  given 
latitude.  This  removes  the  aliasing  effects  of  the  migrating 
tides  with  n=l,3,5,...,  while  leaving  the  effects  caused  by 
the  migrating  tides  with  n  =  2,  4,.... 

[is]  The  second  caveat  is  that  the  different  tidal  modes 
can  cancel  one  another  depending  on  their  phases  and  the 
local  time  of  the  satellite  such  that  the  observed  amplitudes 
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Figure  4.  Altitude-detrended  PLP  plasma  densities  are  shown  as  functions  of  longitude  for  the  night- 
side,  divided  into  2  h  local  time  averages.  The  data  are  averaged  over  3°  wide  longitudinal  bins.  1R1 
model  values  are  shown  for  comparison.  For  most  of  the  night,  there  is  much  weaker  evidence  of  wave¬ 
like  structure,  except  near  dawn  where  a  two-wave  periodicity  can  be  seen.  Deep  density  minima  in  the 
300°-360°  longitude  sector  can  be  seen  in  all  the  plots,  most  clearly  in  the  20:00-22:00  LT  sector. 


Geographic  Longitude 

Figure  5.  Altitude-detrended  PLP  plasma  densities  and 
neutral  density  residuals  are  shown  in  the  20:00-22:00  LT 
sector.  The  neutral  densities  exhibit  wave-4  periodicity. 
Three  of  the  maxima  coincide  with  maxima  in  the  plasma 
density.  The  deep  minimum  in  plasma  density  between 
300°  and  360°  coincides  with  a  similar  deep  minimum  in 
neutral  density. 
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Figure  6.  Altitude-detrended  PLP  plasma  densities  and 
neutral  density  residuals  are  shown  in  the  10:00-12:00  LT 
sector.  The  plasma  densities  exhibit  wave-4  periodicity 
which  is  not  matched  by  the  neutrals. 
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may  exhibit  a  oscillation  at  the  local  time  precession  rate  of 
the  satellite  (or  some  fraction  of  it).  Nevertheless,  from  these 
plots  we  can  glean  information  on  the  nonmigrating  tidal 
structure  on  the  order  of  28  days.  As  the  nonmigrating  tides 
are  highly  variable  with  season,  this  provides  an  attractive 
alternative  to  using  data  from  an  entire  12  h  local  time 
precession  period  ( — 161  days  for  the  GRACE  satellite). 

[19]  Figures  5  and  6  show  the  neutral  density  residuals 
resulting  from  this  analysis,  together  with  the  detrended  PLP 
densities  for  local  time  sectors  20:00-22:00  LT  and  10:00- 
12:00  LT,  respectively. 

4.  Discussion 

[20]  The  average  detrended  plasma  densities  shown  in 
Figures  3  and  4  exhibit  occasional  periodic  structures 
strongly  suggestive  of  tides.  These  periodicities  are  most 
evident  in  Figure  3  (bottom  left),  in  which  four  waves 
appear.  The  same  four-wave  structure  appears  in  Figure  3 
(middle  left),  where  the  IRI  model  also  shows  four  waves. 
Similar  four-wave  periodicity  has  been  shown  in  DMSP 
F13  ion  data  [Ren  et  ah,  2008;  Huang  et  ah,  2010],  CHAMP 
electron  densities  [Liihr  et  al,  2007;  Pedatella  et  al.,  2008] 
and  in  Total  Electron  Content  (TEC)  data  [Lin  et  ah,  2007; 
Scherliess  et  ah,  2008;  Wan  et  ah,  2008;  Pancheva  and 
Mukhtarov,  2010]  and  others.  The  longitudinal  distributions 
of  maxima  and  minima  shown  in  Figure  3  coincide  with 
those  in  the  FUV  observations  [Immel  et  ah,  2006],  in 
agreement  with  predictions  of  the  Global  Scale  Wave  Model 
(GSWM)  [Hagan  et  ah,  2001;  Hagan  and  Forbes,  2002, 
2003].  As  discussed  by  Pedatella  et  ah  [2008]  and  others, 
the  DE3  tide  which  appears  as  a  wave-4  structure  is  dom¬ 
inant  most  of  the  year,  including  the  June  solstice  period. 
It  should  be  noted  that  the  waves  in  the  PLP  data  appear  in 
the  unfdtered  plasma  observations.  Aside  from  the  detrend¬ 
ing  described  above,  the  plasma  data  have  not  been  pro¬ 
cessed  to  extract  wave  modes.  These  tidal  structures  are 
dominant  during  this  period,  thus  the  periodicity  appears  as 
a  macroscopic  feature  of  the  ionosphere,  creating  plasma 
density  variations  up  to  35%  (see  Figure  3,  bottom  left). 
This  sets  our  analysis  apart  from  previous  studies  of  this 
phenomenon. 

[21]  Starting  from  dawn,  shown  in  Figure  3  (top  left),  the 
detrended  plasma  data  show  a  two-wave  structure  that 
evolves  into  a  four-wave  pattern  over  the  following  6  h  of 
local  time.  From  14:00  LT  the  pattern  becomes  less  clear 
with  2  clear  maxima  and  possibly  3  smaller  maxima  evi¬ 
dent.  Note  that  the  data  show  clear  differences  from  the  IRI 
model  in  which  periodicity  is  also  seen  at  08:00-10:00  LT. 
In  general,  the  observed  plasma  density  is  lower  than  pre¬ 
dicted  by  the  model.  In  addition,  the  observed  wave  fea¬ 
tures  deviate  from  the  model  beginning  at  10:00  LT.  As  IRI 
is  an  empirical  model  of  ionosphere,  the  lack  of  wave-4 
periodicity  in  Figure  3  in  the  IRI  traces  shows  that  our  results 
represent  a  new  ionospheric  phenomenon. 

[22]  In  the  nighttime  ionosphere  it  is  considerably  more 
difficult  to  discern  periodic  structure.  Fluctuations  appear  at 
18:00-20:00  LT  in  Figure  4  (top  left)  as  well  as  at  22:00- 
24:00  LT  and  00:00-02:00  LT,  but  they  appear  irregular. 
Shown  on  a  larger  scale,  the  fluctuations  show  some  peri¬ 
odicity,  but  tidal  structures  are  not  obvious.  It  is  possible 


that  these  represent  interference  due  to  multiple  tidal  fre¬ 
quencies,  but  we  have  not  yet  attempted  to  fdter  the  data  to 
determine  if  this  is  the  case.  It  should  also  be  noted  that  the 
C/NOFS  orbital  inclination  of  13°  may  exclude  a  part  of  the 
equatorial  anomalies  where  tidal  features  are  often  detected 
[England  et  ah,  2006;  Liu  and  Watanabe,  2008;  Liihr  et  ah, 
2007].  However,  in  general  the  equatorial  anomalies  are 
weak  during  this  extremely  quiet  solar  minimum  period.  A 
two-wave  pattern  emerges  toward  dawn  (04:00-06:00  LT) 
which  persists  into  the  dayside  (see  Figure  3,  top  left).  This 
is  also  clear  in  the  IRI  model.  For  the  rest  of  the  nightside, 
IRI  results  do  not  display  periodic  structures. 

[23]  Wave-4  periodic  structures  have  been  observed  in 
other  C/NOFS  experiments,  including  the  Ion  Velocity  Meter 
[Araujo -Pradere  et  ah,  2010]  the  Vector  Electric  Field 
Instrument  [Pfaff  et  ah,  2010]  and  density  irregularities  in 
the  PLP  [Dao  et  ah,  2010].  It  is  clear  that  the  ionosphere 
above  400  km  during  solar  minimum  is  dominated  by  tidal 
forces,  generally  thought  to  be  generated  at  low  altitudes 
[Hagan  and  Forbes,  2002,  2003].  As  in  this  study,  all  of 
these  periodic  structures  detected  on  C/NOFS  are  strongest 
on  the  dayside.  This  agrees  with  the  suggestion  by  Immel 
et  ah  [2006]  that  the  upward  propagating  tides  reinforce 
the  regular  E  region  dynamo  which  peaks  at  noon.  However, 
these  results  are  all  obtained  during  a  deep  solar  minimum, 
when  E  region  densities  are  relatively  low.  As  C/NOFS  con¬ 
tinues  to  operate,  we  will  monitor  the  in  situ  densities  for  evi¬ 
dence  of  wave-4  structures  as  solar  activity  increases. 

[24]  The  neutral  density  residuals  obtained  from  the 
GRACE  accelerometer  are  shown  in  Figures  5  and  6  for 
the  dayside  and  nightside,  respectively.  Note  that  the  data 
in  Figures  5  and  6  are  plotted  on  different  scales.  A  com¬ 
parison  of  the  neutrals  with  ions  easily  shows  that  that  the 
tides  are  not  directly  coupled.  In  Figure  5,  the  plasma 
averages  and  neutral  density  residuals  are  plotted  for  the 
20:00-22:00  LT  bin  as  functions  of  longitude.  The  ions  show 
3  maxima  which  coincide  with  maxima  in  the  neutrals,  but 
there  is  little  variation  between  135°  and  270°. 

[25]  In  Figure  6  are  shown  the  ion  and  neutral  densities  for 
the  10:00-12:00  LT  bin.  In  this  case  the  ions  show  wave-4 
periodicity  while  the  neutral  densities  do  not.  The  neutral 
density  residuals  have  been  decomposed  into  longitudinal 
wave  numbers  1-6  (not  shown).  There  is  clear  evidence  of 
wave-4  structure  in  both  the  08:00-10:00  LT  and  10:00- 
12:00  LT  zones.  The  Hough  modes  for  DE3  have  been  com¬ 
pared  with  the  neutral  density  residuals,  and  our  hypothesis 
that  DE3  waves  are  present  has  been  confirmed.  A  full 
description  of  the  neutral  tides  will  be  provided  in  a  separate 
paper. 

[26]  Further  analysis  of  the  neutral  density  residuals  shows 
that  the  results  shown  in  Figure  6  are  due  to  interference 
from  a  series  of  wave  modes.  The  DE3  mode  is  present,  but 
weaker  because  of  diurnal  and  semidiurnal  tides  also  pres¬ 
ent.  In  general  the  neutral  densities  show  wave-4  structure 
most  clearly  on  the  nightside.  The  converse  is  true  of  the 
ion  density  structures.  While  all  tides  have  their  origin  at 
low  altitudes,  their  subsequent  evolution  depends  on  electric 
fields  and  winds  in  the  E  region  which  affect  ion  and  neu¬ 
trals  in  different  ways.  A  causal  link  between  the  wave-4 
structures  and  vertical  plasma  drift  has  been  established 
[Fang  et  ah,  2009].  We  intend  to  explore  this  hypothesis  in 
a  more  comprehensive  study. 
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[27]  One  feature  of  our  analysis  which  is  apparent  in 
both  ions  and  neutrals  is  the  deep  minimum  between  300° 
and  360°  on  the  nightside.  This  feature  corresponds  to  the 
BPDs  initially  identified  [ Huang  et  al.,  2009].  In  the  20:00- 
22:00  LT  sector,  the  average  reduction  in  density  is  85%. 
While  most  obvious  on  the  nightside,  this  minimum  appears 
in  all  of  the  local  time  averages  of  altitude-detrended  ion 
densities  shown  in  Figures  3-6.  It  is  evident  in  the  IRI  values 
shown  in  Figures  3  and  4.  It  is  also  readily  apparent  in 
other  studies  of  ionospheric  tides,  where  it  appears  as  a 
density  decrease  [Fang  et  al.,  2009],  temperature  minimum, 
or  decrease  in  far  ultraviolet  (FUV)  emissions  [Immel  et  al, 
2006].  We  showed  previously  that  the  BPD  is  associated 
with  a  temperature  minimum  in  the  high  F  region  at  DMSP 
altitude  of  840  km  [Huang  et  al.,  2009].  The  simultaneity 
of  this  minimum  in  both  ions  and  neutrals  on  the  night¬ 
side  makes  it  difficult  to  separate  cause  from  effect.  In  situ 
1VM  measurements  of  meridional  drift  show  a  downward- 
directed  velocity  in  the  BPD  (R.  Stoneback,  private  com¬ 
munication,  2010)  which  may  suggest  ion  drag  as  one 
possible  cause  of  the  neutral  density  decrease.  We  believe 
that  ion-neutral  coupling  is  the  cause  of  the  simultaneous 
decrease  in  ion  and  neutral  density  in  this  location,  but 
cannot  postulate  a  physical  mechanism  at  this  time. 


5.  Conclusion 

[28]  We  have  presented  an  analysis  of  altitude-detrended 
plasma  densities  from  the  C/NOFS  satellite  during  June  2008. 
We  used  the  2007  version  of  the  1R1  model  to  remove 
altitude  and  latitude  variations  from  the  observed  densities. 
We  then  averaged  the  results  in  bins  of  local  time  and  lon¬ 
gitude.  This  analysis  revealed  the  presence  of  wave-4  peri¬ 
odic  structures,  strongest  between  08:00  and  14:00  LT.  On 
the  nightside,  the  periodicity  is  much  weaker,  but  there  is 
a  persistent  deep  minimum  in  plasma  density  at  approxi¬ 
mately  330°. 

[29]  We  analyzed  the  neutral  density  variations  from  the 
GRACE  satellite  during  the  same  period,  removing  altitude 
variations  and  restricting  the  data  to  13°  in  latitude  from  the 
equator,  similar  to  the  C/NOFS  orbit.  We  averaged  the 
density  residuals  in  similar  local  time  and  longitude  bins  as 
for  the  ion  densities.  The  results  exhibit  wave-4  periodicity, 
strongest  on  the  nightside.  This  is  the  converse  of  the 
ionospheric  periodic  structures  which  are  strongest  on  the 
dayside. 

[30]  We  conclude  that  periodic  structures  in  both  ion  and 
neutral  densities  are  probably  generated  by  tides  at  low 
altitudes.  However,  their  subsequent  evolution  from  their 
common  source  results  in  a  decoupling  of  the  structures  at 
F  layer  altitudes.  One  striking  exception  is  the  large  decrease 
in  density  which  is  seen  in  both  ions  and  neutrals  at  330°  on 
the  nightside.  We  believe  this  to  be  the  result  of  ion-neutral 
coupling,  but  a  mechanism  to  create  this  density  decrease  is 
beyond  the  scope  of  this  report. 
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